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ABSTRACT 

Diffusive cosmic-ray transport in nonuniform large-scale magnetic fields in 
the presence of boundaries is considered. Reflecting and absorbing boundary 
conditions are derived for a modified telegraph equation with a convective term. 
Analytical and numerical solutions of illustrative boundary problems are pre¬ 
sented. The applicability and accuracy of the telegraph approximation for focused 
cosmic-ray transport in the presence of boundaries are discussed, and potential 
applications to modeling cosmic-ray transport are noted. 

Subject headings: cosmic rays — diffusion — magnetic helds — scattering — 

Sun; particle emission 


1. Introduction 

When energetic cosmic-ray particles propagate in the cosmos, they often interact with 
turbulent magnetic helds. The evolution of a particle distribution is governed by the Fokker- 
Planck equation (e.g., Schlickeiser 2011, and references therein). When pitch-angle scattering 
is strong enough to ensure that the scale of density variation is signihcantly greater than the 
particle mean free path, the evolution can be approximated as a simpler diffusive process 
(Hasselmann & Wibberenz 1970; Schlickeiser & Shalchi 2008). 
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The diffusion approximation has been employed for studying the acceleration and prop¬ 
agation of energetic particles in various astrophysical situations—from the atmosphere of the 
Sun and interplanetary space (e.g., Bieber et ah 1987; Le Roux & Webb 2009; Artmann et 
ah 2011) to the interstellar medium (e.g., Schlickeiser, 2009; Litvinenko & Schlickeiser 2011; 
Schlickeiser et ah 2011). Observations reveal large-scale spatial variations of the magnetic 
held in all those locations (e.g., Sofue et al. 1986; Sandroos & Vainio 2007; Droge et ah 
2010). The diffusion model has been developed to incorporate the coherent particle stream¬ 
ing due to the effect of adiabatic focusing in a nonuniform background magnetic held (Earl 
1981; Beeck & Wibberenz 1986; see also Litvinenko 2012a, 2012b; He & Schlickeiser 2014, 
and references therein). 

Propagation of solar energetic particles in interplanetary magnetic helds remains a sub¬ 
ject of intense research activity (e.g., Zhang et al. 2009; Droge et al. 2010; Qin et al. 
2013; Laitinen et al. 2013; Wang & Qin 2015). The theoretical studies are motivated by 
the new data from multi-point spacecraft observations, which allow new insights into the 
physics of the particle transport (e.g., Dresing et al. 2014; Droge et al. 2014; Lario et al. 
2014). Since theoretical modeling is usually based on numerical solutions, simple analytical 
approximations can complement and guide the simulations. 

A general shortcoming of the dihusion approximation is that the dihusion equation 
implies an inhnite speed of signal propagation, whereas particle speeds are hnite, of course. A 
more accurate description may be provided by the telegraph equation (Goldstein 1951). This 
is an equation of hyperbolic type, and its solution at long times asymptotically approaches 
the solution of the diffusion equation (Davies 1954). The telegraph equation for Brownian 
motion (Brinkman 1956; Sack 1956) had been shown to be substantially more accurate than 
the diffusion equation (Hemmer 1961). The derivation of a generalized telegraph equation 
and its applications for modeling cosmic-ray transport had been repeatedly considered both 
in the limit of a uniform background magnetic held (e.g., Fisk & Axford 1969; Earl 1974, 
1992; Gombosi et al. 1993; Schwadron & Gombosi 1994) and in the more realistic case of 
a spatially varying held (e.g., Earl 1976; Pauls & Burger 1994; Litvinenko & Noble 2013; 
Litvinenko & Schlickeiser 2013; Maikov & Sagdeev 2015). 

Boundaries typically need to be considered in cosmic-ray transport problems in both 
interplanetary and interstellar plasmas (e.g., Schlickeiser 2009; Artmann et al. 2011). The 
boundary conditions for both the dihusion equation and the standard telegraph equation 
are well known (Masoliver et al. 1992, 1993). We are not aware of a published derivation of 
the corresponding boundary conditions for focused particle transport, described by a modi- 
hed telegraph equation. In this paper, we develop the telegraph approximation for focused 
cosmic-ray transport in the presence of boundaries. The new analytical and numerical results 
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complement those for cosmic-ray transport in the absence of boundaries, which we presented 
in our previous study of the telegraph approximation (Effenberger & Litvinenko 2014). 

Following a brief description of the model (Section 2), we derive the reflecting and 
absorbing boundary conditions (Section 3), illustrate the use of the Laplace transform and 
Fourier series for obtaining analytical solutions (Section 4), and compare them with numerical 
solutions of the Fokker-Planck equation (Section 5). 


2. The telegraph approximation for focused cosmic-ray transport 


Spatial non-uniformity of the mean magnetic held leads to the adiabatic focusing that 
results in coherent streaming of cosmic-ray particles along the mean held (Roelof 1969; Kun- 
stmann 1979). Earl (1976) derived a modihed telegraph equation for the focused particle 
transport in a spatially varying magnetic held. The equation, however, described the co¬ 
efficient of an eigenfunction expansion rather than the particle density that is the physical 
quantity of interest. The problem was recently reexamined, and the telegraph approximation 
was derived for the particle density in a spatially varying magnetic held in a weak focusing 
limit (Litvinenko & Schlickeiser 2013). The telegraph approximation has also been obtained 
in a complementary case of an arbitrary constant focusing strength and isotropic pitch-angle 
scattering (Litvinenko & Noble 2013), using an exact expression for the variance of a particle 
distribution, obtained by Shalchi (2011). 


A possible physical context for the telegraph approximation is provided by the Fokker- 
Planck equation for focused transport, which describes the evolution of the distribution 
function of energetic particles: 


dfo 

dt 


+ /iU 


dfo 

dz 


+ 


V 




dfo 

(9/i 


A 

(9/i 



( 1 ) 


(Roelof 1969; Earl 1981). Here fo{z, is the distribution function (gyrotropic phase- 

space density), t is time, jj, is the cosine of the particle pitch angle, v is the particle speed, 
2 ; is the distance along the mean magnetic held Bq, L = —Bo/^dBo/dz) is the adiabatic 
focusing length, and is the Fokker-Planck coefficient for pitch-angle scattering. 


To illustrate the application of the approximation to a model transport problem, we 
consider the simplest physically plausible model. Specihcally, we assume a constant focus¬ 
ing length L, and we neglect momentum dihusion, adiabatic cooling, and advection with 
a background how (say, the solar wind). The simplifying assumptions are discussed, for 
example, by Artmann et al. (2011) and Ehenberger & Litvinenko (2014) in the context of 
solar energetic particle transport in interplanetary space. 



4 


The distribution function can be expressed as the sum of the isotropic density Fq and 
an anisotropic component Qq: 


fo{z, /i, t) = Fo{z, t) + Qq^Z, fi, t), 


( 2 ) 


where 

1 

Fo{z,t) = -J fodfx. (3) 

Assuming that go Fq, an approximate expression for go can be found and substituted into 
Equation Q, integrated with respect to /r. Depending on the accuracy of the expression for 
go, the result is either the usual diffusion approximation or the telegraph approximation for 
focused transport (e.g., Earl 1976; Litvinenko & Schlickeiser 2013). Here we investigate the 
modihed telegraph equation for the isotropic density: 


dFo 


+ T- 


d^Fo 


dt ' ' dF ’ 

where is the coherent speed (e.g., Earl 1981), Ky is the parallel diffusion coefficient, 
^ = Xo/L is the focusing strength, and Aq is the scattering mean free path in the absence 
of focusing (.^ = 0). Equation (|^ is written in dimensionless form by measuring distances 
in units of Aq, speed in units of the constant particle speed v, and time in units of Xo/v. 
Although we formally recover the diffusion approximation by setting r = 0, in practice r is 
of order unity in a physically relevant parameter range (Litvinenko & Schlickeiser 2013). 


d^Fo 


+ 


dFo 


( 4 ) 


The coefficients k\\ and r generally depend on the focusing strength In the weak focus¬ 
ing limit 1, the transport coefficients are given by Equations (9) and (14) in Litvinenko 
& Schlickeiser (2013) for an arbitrary pitch-angle scattering coefficient For isotropic 

scattering, the telegraph equation has been derived for an arbitrary focusing strength 
and K\\ and r are given by Equations (9) and (24) in Litvinenko & Noble (2013). The 
complementary expressions agree in the limit of weak focusing and isotropic scattering. 

It is worth noting that the derivation of the telegraph equation by Earl (1973) and 
Litvinenko & Schlickeiser (2013) involves truncating an inhnite system, equivalent to the 
original Fokker-Planck equation. A well-known weakness of that approach (Gombosi et ah 
1993; Schwadron & Gombosi 1994) has been recently reiterated by Maikov & Sagdeev (2015). 
Truncating the system of harmonics, derived from the Fokker-Planck equation, leads to an 
error in the coefficient r—and more generally, to an equation which is not correctly ordered. 
The correct way of deriving the telegraph equation relies on an asymptotic expansion that 
yields the diffusion equation in the lowest order and the telegraph equation in the next order. 
Gombosi et ah (1993) described the method in detail for the case of a uniform background 
magnetic held and isotropic scattering. Of course the correct mathematical procedure is gen¬ 
erally preferable. In practice, however, the improvement in accuracy may be quite modest. 
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For instance, Gombosi et al. nsed the parameter r to calcnlate the signal propagation speed 
in the telegraph eqnation by both methods. The correct derivation gave n-^5/11, where v 
is the particle speed, compared with when simple trnncation was nsed—abont 17% 

difference. Effenberger and Litvinenko (2014) showed that the telegraph eqnation for fo- 
cnsed transport, obtained by the simpler method, yields a solntion that agrees well with the 
nnmerical solntion of the original Fokker-Planck eqnation on an infinite interval. It is also 
worth stressing that the telegraph approximation was shown to reprodnce the evolntion of 
the particle density profile mnch more accnrately than the diffnsion approximation (r = 0), 
especially when the focnsing is strong (Litvinenko & Noble 2013; Effenberger & Litvinenko 
2014). 

It may be nsefnl to rewrite Eqnation ([^ in terms of the linear density 

F(z,t) = exp(^z)Fo(z,t), (5) 


which is the nnmber of particles per line of force per nnit distance parallel to the mean 
magnetic held. Clearly, the linear density satishes 


OF 


d‘^F 


d‘^F 


OF 

The descriptions in terms of Fq and F are mathematically eqnivalent, and the choice of Fq or 
F is a matter of convenience (Earl 1981). For instance, particle conservation in the absence 
of sources and sinks is conveniently expressed as 


, OF 


( 6 ) 


N(t) = 2 


Fdz = const. 


( 7 ) 


Note for clarity that Litvinenko & Schlickeiser (2013) did not work with the linear density 
and they used the notation F{z,t) instead of Fo(^, t). 


3. Boundary conditions 

Masoliver et al. (1992, 1993) obtained boundary conditions for the standard telegraph 
equation in the presence of rehecting or absorbing boundaries. We generalize the arguments 
in Masoliver et al. (1992, 1993) and derive the boundary conditions for a modihed telegraph 
equation with a convective term that can describe the focused cosmic-ray transport. 

Decomposing the particle distribution function into two components corresponding to 
particles moving to the right, a{z,t), and to the left, b{z,t), we have the coupled equations 
that reduce to Equation (1) in Masoliver et al. (1992) in the limit ^ = 0: 

dttt = + —{b — a) + 


( 8 ) 
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dtb = 



(9) 


The linear density 


F{z, t) = a + b 


( 10 ) 


satisfies the focused telegraph equation (j^. Physically, w = is a signal propagation 

speed in the telegraph equation, which does not depend on the first-order terms in the 
equation. The first term on the right in Equation (|^ describes convective transport, the 
second term describes the particle change of direction, and the third term describes the 
focusing effect of a nonuniform background magnetic field. The interpretation is similar for 
Equation (|^. 


Consider first a region zi < z < Z 2 with absorbing boundaries at z = zi and z = Z 2 , so 

that 

a{zi,t)=0, b{z 2 ,t) = 0. (11) 

Subtracting Equation ([^ from Equation ([^ and using Equation ( [II| yields the boundary 
conditions 

± ^{d,F -^F) = F + TdtF, (12) 


where the plus (minus) sign corresponds to the left (right) boundary at z = Zi {z = Z 2 ). The 
absorbing boundary conditions for the density Fq follow from Equations ([^ and (12): 


± = Fq + rdtFo. 


(13) 


If we formally set r = 0, the telegraph approximation simplifies to the focused diffusion 
model, termed pseudo-diffusion by Earl (1981), and the absorbing boundary conditions above 
simplify to the familiar condition Tq = F = 0 at an absorbing boundary. 

Now suppose that reflecting boundaries are present at z = Zi and z = Z 2 , so that 


a{zi,t) = b{zi, t), a{z 2 , t) = b{z 2 , t). 


(14) 


In this case, subtracting Equation ([^ from Equation 
the boundary conditions 

d,F - = 0 


and 


d.Fc) = 0 


and using Equation (14) leads to 

(15) 

(16) 


at z = zi and z = Z 2 , which are the same conditions as those for both the diffusion equation 
and the telegraph equation considered by Masoliver et ah (1993). In the case of a single 
reflecting boundary. Equation (15) immediately follows on integrating Equation (§ and 
using the particle conservation given by Equation ([^. 
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Finally, we note that Marshak (1947) advocated the use of an escape condition for the 
particle flux 




ligodfi 


(17) 


'-1 


as a more accurate alternative to the conventional absorbing boundary condition. Marshak’s 
condition is given by 


at the left boundary and 






fivFodfi 


'-1 


fivFQdfi 


(18) 


(19) 


at the right boundary. For example, in the weak focusing limit of the telegraph approxima¬ 
tion, the flux is given by Equation (13) in Litvinenko & Schlickeiser (2013), which yields 


- Kiid^Fo TKiidtzFo = ±-Fo 


( 20 ) 


in our dimensionless units, where the minus (plus) sign corresponds to the left (right) bound¬ 
ary at z = zi {z = Z 2 ). In the limit r = 0, the condition reduces to the corresponding bound¬ 
ary condition of the diffusion model (Weinberg & Wigner 1958). In practice, the presence 


of a second mixed derivative dtzFo in Equation (20) makes the boundary condition difficult 
to use. 


4. Analytical solutions 
4.1. Infinite interval 

Analytical solutions of the focused diffusion model are commonly used to model the 
transport of solar energetic particles in interplanetary space (e.g., Artmann et al. 2011 
and references therein). Since the telegraph approximation should give a more accurate 
description of the cosmic-ray transport, analytical solutions of the telegraph approximation 
could be used to quantify the accuracy of the diffusion approximation or to validate the 
accuracy of a numerical solution of the Fokker-Planck equation for an evolving particle 
distribution function. Here we illustrate how initial and boundary value problems for the 
modified telegraph equation can be solved by Laplace transform for infinite and semi-infinite 
intervals and by Fourier series for a finite interval. 

For simplicity we assume a symmetric point source at z = zq, so that a(z, 0) = b{z, 0), 
and the initial conditions are given by 

Fo(2.0) = <S(s - Jo), d,F„{z,0) = 0, 


( 21 ) 



where the second condition follows from Equations ([^ and (|^. The telegraph equation is 
linear, and so we omit the normalization constant | exp(—iii -^ 0 (^ 50 ) for brevity. 

Consider first the initial value problem for an infinite interval. Although in this case 
the solution can be obtained by a change of variables, which reduces Equation Q to the 
standard telegraph equation (e.g., Kevorkian 2000), it is instructive to solve the problem 
directly using the Laplace transform. The transform of Equation Q reads 

Fq + ^Fq -s(l + ts)Fq = -(1 + ts)5{z — zq), (22) 

K|| Ky 


with the Laplace transform Fo{z,s) = L[Fo(z,f)] and the prime denotes differentiation with 
respect to 2 ;. Solving the ordinary differential equation yields the transform 


Fo{z,s) = 


1 + rs 
2K\\r] 


exp ( -|(z - zo) - r]\z - Zq 


where 




s(l + TS) 


K\\ 


Consequently the solution of the initial value problem is given by 

Fo{z, t) = Go + rdtGo, 

where 


Go{z, t) — L 


-1 


i-^{z - Zo) - ri\z - Zo 


.2«||h 


(23) 


(24) 


(25) 


(26) 


is the fundamental solution of the modihed telegraph equation. Evaluating the inverse 
transform yields 


Go{z,t) = ;^-4=exp (-^{z-zo) - Io(u)F (t- \ l 'Iz - zo 


2^jqf 



ay 


(27) 


Here, H is the Heaviside step function, Jq is a modified Bessel function, and its argument is 

1 


u 


" /C W “ ^0)^ 




H 


T 


(28) 


Note that we assume C^«^yT < 1. The inequality was shown to be valid both in the case 
of weak focusing and anisotropic scattering (Litvinenko & Schlickeiser 2013) and in the 
case of strong focusing and isotropic scattering (Litvinenko & Noble 2013). For instance. 












(1 — = 1 for isotropic scattering. Substituting Equation (27) into Equation (25) 

leads to a solution in terms of the modified Bessel functions Jq and Ii. 




/o(m) + (1 

5 


t h{u) 
2t u 


H t- z 



Zo\ 



- [z - Za) \ + S IJ+ {z - z„) 


(29) 


which slightly generalizes Equations (18) and (19) in Effenberger & Litvinenko (2014) who 
compared the solution with both the prediction of the diffusion model and a numerical 
solution of the original Fokker-Planck equation. 


4.2. Semi-infinite interval: a reflecting boundary 


Now consider the initial value problem, specifled by Equations (16) and (21), on a 


semi-inhnite interval z > 0. Physically, a reflecting inner boundary condition at = 0 may 
correspond to the transport of solar energetic particles, accelerated close to the solar surface. 
The reflection of particles traveling towards the sun is caused by a magnetic bottle effect of 
the strongly converging magnetic held. In this case the solution of the transformed Equation 


(22), which satishes dzFo{0,s) = 0, is given by 


’(+ exp{-ri\z - Jol) + (>)- tq exp(-))(z + Jo)) 


(30) 


As an interesting aside, note that the mean age T of particles at a given location z can 
be elegantly expressed in terms of the Laplace transform Fq: 


T = 


I^tFodt 

J-F,dt 


= -dsFo{z,0). 


(31) 


For example, suppose that a particle source is located very close to the solar surface, so that 
we may take Zq = 0. We calculate 

1 Z 


which agrees with the advection-dominated limit of Equation (4) in Jokipii (1976). 


( 32 ) 
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The solution of the initial value problem is again given by Equation (25), but now the 
fundamental solution is as follows: 


Go{z, t) — L 


-1 




exp {-r]\z -zo\) + I exp {-r]{z + zq)) 
2t] + ^ 


(33) 

For simplicity, consider again a particle source located near the boundary. Setting Zq = 0 
and keeping in mind that z > 0, we have 


Go{z, t) — L 


-1 




(34) 


The inverse Laplace transform can be expressed in terms of a table transform (Erdelyi et ah 
1954) by noticing that 

r 


r /2 = 


K|| 




where 


2 ^ 1 - ? 'tr 

4t2 ■ 


(35) 


(36) 


The resulting fundamental solution is as follows: 




exp 1 1 +« 




exp 


— — w'^ Ii{aw)dw 


(37) 


The result can be verihed in the limiting case of a uniform magnetic held (.^ = 0) when 
the solution simplihes to 


G„(..i) = ^exp(-^|/„ 




where 


Consequently, 


1 z'^ 

2yr2 Kyr' 




Io{Uo) + 


t Iljup) 
2t uq 


Hit- z 





(38) 


(39) 


(40) 
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In other words, the solution on the semi-inhnite interval z > 0 with a source at the edge 
of the interval in the absence of focusing is simply double the solution of the corresponding 
initial value problem for an infinite interval. The result immediately follows on applying the 
method of images to the problem. 

Finally, for a uniform magnetic field, the solution on a semi-infinite interval with an 
absorbing boundary is given by Equation (19) in Masoliver et al. (1992). 


4.3. Finite interval: reflecting boundaries 


The relatively simple form of the reflecting boundary conditions makes it possible to 
express the solution of the telegraph equation on a finite interval in terms of a Fourier series. 
The series solution is particularly convenient for long times, when the first few terms of 
the series accurately approximate the solution. Suppose that two reflecting boundaries are 
present at 2:1 = 0 and Z 2 = I, Fo{z, 0) is given for 0 < 2 ; < Z, and dtFo{z, 0) = 0 for simplicity. 
Straightforward application of the method of separation of variables to Equations (|^ and 


(16) leads to 


Fo{z,t) = Co + ^CnZn{z)Tn{t), 

n=l 


where 


ZJz) = 


z) = exp 


' 2 


cos 


/Tmz\ 

) + 


V I 


El /"irnz 

-sm —— 

27rn V I 


(4) 

cosh OJnt -\ -sinh OJnt 

V 2ry 

20JnT 


UJn = 


1 

47-2 


K|l / /Trn\‘^ 


11/2 


The initial conditions yield 


(41) 

(42) 

(43) 

(44) 


Co = 




[expiEl) - 1 ] 


exp (^z) Fq{z, 0)dz, 


(45) 


Cn 



/"Knz\ El /7mz\ 

cos —— -I -sin —— 

\ I ) 2Tm \ I ) 


Fq{z, Q)dz. 


(46) 


In the limit ^ = 0, the solution for the initial profile Fq[z, 0) = S(z — zo) agrees with Equation 
(19) in Masoliver et al. (1993). 
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4.4. Finite interval: absorbing boundaries 


The appearance of two partial derivatives in Eqnation (13) complicates the solntion of 


the telegraph eqnation with absorbing bonndaries. In the case = 0 and -Fo(z, 0) = 6{z — Zo), 
Eqnation (25) in Masoliver et ah (1992) gives an exact solntion in terms of an inhnite 
series of Bessel fnnctions. Since nnmerical solntions might not captnre important qnalitative 
featnres of the solntion, it is natnral to seek an approximate analytical solntion. Here, we 
illnstrate an integral approximation method, similar to the heat-balance approximation in 
heat condnction (e.g.. Crank 1984; Hill & Dewynne 1987). 

Consider a hnite region with absorbing bonndaries at 2:1 = 0 and Z 2 = I- The bonndary 
valne problem for Eqnation (|^ might serve as the basis for a simple model for the escape 
of galactic cosmic rays away from the galactic plane (Schlickeiser 2009). The idea of the 
approximation is that, instead of solving Eqnation (|^ exactly, we integrate it with respect 
to 2 ; from 0 to I and seek a solntion of a specihed fnnctional form that satishes the bonndary 
conditions. As a simple illnstration, consider 

Fo{z,t) ^ fi{z)f 2 {t) (47) 

and reqnire that it satishes both the absorbing bonndary conditions, given by Eqnation ( |I^ , 
and an ordinary differential eqnation for / 2 (t), obtained by integrating Eqnation (|^ over 2 ;. 
On choosing a parabolic density prohle 

fi{z) = 1 + kz{l - z), (48) 

it is straightforward to verify that 

rh + fy + , , L,j 2 = 0, (49) 


1 + kP/Q 

and so a possible solntion is given by 

/ 2 (f) = const exp(Af), 


where the constants k and A follow from Eqnation (49) and the bonndary conditions: 

1 -|- At 


k = 


It follows that 






K|| 1 



r 1 -|- kP/Q 



(50) 


(51) 


(52) 


(53) 


and so the approximate solntion describes the evolntion of an initially broad (parabolic) 
particle density prohle. 
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5. Stochastic simulations of the Fokker—Planck equation 
5.1. Numerical method 

To quantitatively assess the accuracy of the telegraph approximation in the presence of 
boundaries, we use numerical solutions of the corresponding Fokker-Planck equation. The 
numerical approach is based on solutions to an equivalent system of stochastic differential 
equations (SDEs), similar to the method employed in Litvinenko & Noble (2013) and Effen- 
berger & Litvinenko (2014). In the following, we use the Milstein approximation scheme, 
given by Equations (29) and (30) in Effenberger & Litvinenko (2014), and we refer the reader 
to this paper for more details on the numerics (see also Kopp et ah 2012, and references 
therein for more details on the SDE method). 

In the present study, we have to take special care of the boundary conditions required for 
the Fokker-Planck equation. In practice, this means that the trajectories of pseudo-particles 
are integrated according to their stochastic evolution equations, until they cross a boundary 
(say aX z = Q oi z = 1). At this point, either the particle speed is reversed (and so its pitch 
angle changes its sign at a reffecting boundary) or the particle is discarded for the rest of 
the simulation at an absorbing boundary. The final distribution function is constructed in 
the usual way, by applying an appropriate binning procedure to the particle positions, and 
normalizing the result for comparison with the analytic predictions. 

We consider the case of isotropic scattering, for which the Fokker-Planck pitch-angle 
scattering coefficient is given by 

= Do(l-/i^)- (54) 

In this case, the coefficients in the telegraph equation (|^ are given by 

coth ^ 1 

-11 = — 

tanh,^ 


(Litvinenko & Noble 2013; He & Schlickeiser 2014, and references therein). We have Ky ~ 1/3 
and r ~ 1 in the weak focusing limit < 1. 


(55) 

(56) 


5.2. Semi-infinite interval: a refiecting bonndary 

We first consider the positive semi-infinite region with a refiecting boundary at 2 ; = 0. 
Particles are injected at time to = 0 and position = 1 with isotropic pitch angle. Masoliver 
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Fig. 1.— Isotropic density F{z,ti) at four different times = 1 (black, ’x’), t 2 = 2 (red, 
’+’)! ^3 = 5 (blue, ’o’) and = 10 (green, ’o’) for the semi-infinite domain with reflecting 
boundary. The solid lines show the solution of the telegraph equation, given by Equation 
(57). The symbols show the numerical results, obtained from 10® particles starting at Zq = 
1 . 0 , and averaging without regard to the pitch angle of the particles. 


et ah (1993) gave a solution for the telegraph equation for a case without focusing based 
on the method of images (their Equations (14) and (16)). Their solution follows from our 
Equation (29) in the limit 0 when we apply the method of images: 

r^\z,t) = Fo{z,t\zo) + Fo{z,t\-zo), (57) 


where Fo(z,t|zo) denotes the solution with injection at Zq. Note for clarity that in Masoliver 
et ah (1993) the signal propagation speed is denoted as c and their T is equal to our r. The 
method of images is not applicable in the case ^ 7 ^ 0 . 

We concentrate on the case .^ = 0 to describe the effect of the boundaries on the par¬ 
ticle distribution. We compare the analytical solution for the semi-inhnite domain with 
simulation results for the Fokker-Planck equation as described above. Figure [T] shows the 
evolution of the particle distribution at four different times, calculated from an injection of 
10® pseudo-particles. The comparison shows a good agreement after a few scattering times, 
as was the case for the solution on an inhnite interval (Effenberger & Litvinenko 2014). 
The disagreement between the numerical and approximate analytical solutions, however, is 
signihcant at earlier times. The overall amplitude of the distribution appears to be under¬ 
estimated by the telegraph solution. Note that, while the telegraph solution conserves the 
particle number, the ^-functional contributions are omitted in Figure [T] Furthermore, the 
slower signal propagation speed, when compared to the particle speed, results in positions of 
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the discontinuities, which underestimate the fast particle spread beyond the limits imposed 
by the signal propagation speed, especially at early times. At times t = 5 and t = 10 the 
agreement is signihcantly better, and the discontinuities have almost disappeared. 



Fig. 2.— Time prohle at ^ = 1 for the semi-inhnite interval with a reflecting left boundary. 
Particles are injected at zq = 0.5 to illustrate the reflection effect and the resulting discon¬ 
tinuities in the telegraph solution. The solid black line gives the solution of Equation (57), 
while the dashed line gives the solution for an inhnite interval (without a reflecting bound¬ 
ary at the origin). The symbols connected with a red line show the numerical SDE solution 
obtained from the injection of 10® particles. 


In addition to the spatial density prohles, time prohles of particle intensities can provide 
an important tool for analyzing solar energetic particle data. Figure illustrates the effect 
of a reflecting boundary at z = 0, say due to the magnetic mirroring of energetic particles 
accelerated in the solar corona. To emphasize the effect of the reflecting boundary for il¬ 
lustrative purposes, we assume that the particle injection occurs at zq = 0.5 and plot the 
time prohle at z = 1. A jump in the particle intensity caused by the arrival of rehected 
particles is clearly visible. The early arrival of non-scattered particles at dimensionless time 
t = {z — Zq)/v = 0.5 causes the discrepancy between the numerical Fokker-Planck solution 
and the telegraph approximation, which almost disappears when the diffusing rehected par¬ 
ticles start arriving at time t = [z t- zq)/w k, 2.6 and decreases even further as the transport 
becomes more dihusive for later times. 
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Fig. 3.— Similar to Figure [T| but now with an absorbing left boundary condition. 

5.3. Semi-infinite interval: an absorbing boundary 


We consider a setting similar to the one in Section 5.2, with initial particle injection again 
at Zq = 1, but now for a left-hand side absorbing boundary condition on the semi-inhnite 
interval. Figure [^compares the analytic result for the case without focusing (Equation (19) 
in Masoliver et ah 1992) with the numerical Fokker-Planck results for 10® particles. For 
the purpose of comparison, the solution of Masoliver needs to be rescaled with z = cy and 
normalized with the signal speed, due to the different nondimensionalization. 


As in the previous case with reflecting boundary conditions, we hnd that the agreement 
of the analytical and numerical solutions improves as later times. A noteworthy feature is 
the non-vanishing boundary values at an absorbing boundary, in contrast to the boundary 
condition F = Fq = 0 in the diffusion approximation (see also Section 5.5 below). 


5.4. Finite interval: reflecting boundaries 


As the next example, we consider a hnite interval with reflecting boundaries. We eval¬ 


uate equation (41) with N = 1000 terms in the Fourier series on a hnite domain of length 
I = 10 with an isotropic injection of particles at Zq = 2.5. Figure shows a comparison 
between this result and a stochastic simulation of 10® particles giving a solution to the cor¬ 
responding Fokker-Planck Equation ([^. We consider the cases of no focusing (^ = 0) and 
strong focusing (^ = 1.5). 


At relatively early times, say ti = 5, the Fourier series solution shows limited applicabil- 
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Fig. 4.— No focusing (.^ = 0, left panel) and strong focusing = 1.5, right panel) results 
for the isotropic linear density F{z,ti) on a hnite domain of length I = 10 with injection 
at Zq = 2.5 and reflecting boundaries. The four different times displayed are ti = 5 (black, 
’x’), t 2 = 10 (red, ’+’), = 20 (blue, ’o’) and = 50 (green, ’o’). The solid lines show the 

solution of the telegraph equation, given by Equation (41). The symbols show the numerical 
results, obtained from 10® particles starting at zq. In the left panel, the dashed line indicates 
the constant steady-state solution. 


ity near the discontinuities in the solution of the telegraph equation. Similar to the previous 
cases of inhnite and semi-inhnite domain, the Fokker-Planck solution exhibits a relatively 
smooth prohle at the discontinuities, and at later stages the agreement is very good, both for 
weak and strong focusing. For long times, the solution approaches the steady state solution. 


given by the term cq in Equation (41): 


Fq(z, t —^ Oo) —^ Cq — 

which reduces to Fq = N/21 in the limit ^ = 0. 




2(exp{^l) - 1)’ 


(58) 


5.5. Finite interval: absorbing boundaries 

Finally, we consider a hnite domain of length I = 10 with two absorbing boundaries. We 
compare the approximate parabolic solution of section 4A to SDE solutions of the Fokker- 
Planck equation, initialized either directly with the parabolic prohle of equation (48) or with 
a point injection at Zq = 5. The values for k and A are determined from Equations (51) and 
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Fig. 5.— Comparison of the parabolic solution for the finite interval with absorbing bound¬ 


aries (Equation (47), solid lines) to a SDE solution (symbols) of the Fokker-Planck equation 
with 10® particles initialized with the parabolic profile at f = 0 (solid black line) with a 
rejection sampling method (left panel) or with a point injection ai zo = 5 (right panel). The 
three different times displayed are = 1 (red, ’-|-’), t 2 = 10 (blue, ’o’) and = 20 (green, 
’o’). 


(53) as 


A = -0.0304 (59) 

k = 0.1679. (60) 

The second possible value for A would result in a negative value for k giving an inverted 
parabolic prohle. 

Figure]^ shows the two cases (parabolic initial prohle, left panel; delta-functional injec¬ 
tion, right panel). For the left panel plot, the solutions have been normalized to agree for 
the initial prohle. For the right panel, the solution is normalized to agree approximately at 
time fs = 20. It is clear that the parabolic solution captures both the general structure of 
the Fokker-Planck solution and its approximate exponential decay behavior quite well. The 
detailed time evolution and spatial dependence, however, are more complicated, and for the 
delta-functional injection, the parabolic prohle only develops at a later stage. As mentioned 
before, in the context of the semi-inhnite interval, the telegraph approximation allows for 
non-vanishing boundary values, in contrast to the boundary condition F = Fq = 0 in the 
dihusion approximation. This appears to yield good agreement of the telegraph solution 
with the numerical Fokker-Planck results in the vicinity of the absorbing boundaries. 
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We could not confirm numerically the presence of discontinuities in the analytical so¬ 
lution of Masoliver et ah (1992) for the initial isotropic pitch-angle distribution of particles 
(see Dunkel et ah 2007, for discussion of this feature of the telegraph approximation, as 
well as for an alternative approach to diffusive transport modeling). An initial condition 
corresponding to two oppositely-traveling particle beams might give a better agreement of 
the numerical Fokker-Planck solution and the telegraph approximation. 


6. Discussion 

In this paper we systematically developed the telegraph approximation for particle trans¬ 
port in the presence of boundaries, taking into account the adiabatic focusing effect in a non- 
uniform mean magnetic field, which leads to coherent particle streaming along the held. We 
derived rehecting and absorbing boundary conditions for a modihed telegraph equation with 
a convective term, and we presented analytical solutions of illustrative boundary problems, 
which might be relevant for modeling diffusive cosmic-ray transport in nonuniform large-scale 
magnetic helds in the presence of boundaries. We also demonstrated the accuracy of the 
telegraph approximation for focused transport in the presence of boundaries by comparing 
the analytical solutions of the telegraph approximation with the numerical solutions of the 
original Fokker-Planck equation. The numerical results complement those for an inhnite 
interval (Litvinenko & Noble 2013; Effenberger & Litvinenko 2014). 

The key point in assessing the practical usefulness of the telegraph model is that cosmic- 
ray transport in turbulent magnetic helds is dihusive on sufficiently long time scales. After a 
few scattering times, the dihusion and telegraph models give very similar predictions for an 
evolving particle distribution. The telegraph approximation is an attempt to more accurately 
describe the particle evolution on shorter time scales. Direct comparison of the predictions of 
the dihusion and telegraph models with the numerical solution of the Fokker-Planck equation 
for focused particle transport clearly shows that the telegraph model reproduces the shape 
of an evolving density pulse much better than the dihusion model, especially when focusing 
is strong, even for times signihcantly exceeding the scattering time (see, for instance. Figures 
2 and 3 in Ehenberger & Litvinenko 2014). 

In the present study, we compared the Fokker-Planck and telegraph results for diherent 
boundary value problems and conhrmed the validity of the telegraph approximation after 
just a few scattering times. Hence the telegraph approximation can provide an improved 
description over the dihusion approximation at these early times. Values for the parallel 
mean free path of about 0.1 to 0.3 AU, found in solar energetic particle studies (e.g. Droge 
et al. 2014), suggest that the time scale of a few scattering times is relevant for interplanetary 
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particle transport, for example in the problem of predicting the impact of large solar events 
at 1 AU (e.g. Shea & Smart 2012). 

A limitation of the telegraph approximation is that solving an initial value problem 
requires the knowledge of the first derivative of the density FQ{z,t) with respect to time at 
t = 0, which can be determined only if the full distribution function /q^z, is known. 

However, in concrete applications OtF^lz, 0) may be known or estimated either observation- 
ally or on theoretical grounds. For instance, in the simplest case of an isotropic initial 
distribution, the derivative vanishes (see equations 4 and 5 in Litvinenko and Schlickeiser 
2014). More generally, any knowledge of the initial distribution would yield a more accurate 
description of the evolving particle distribution in the telegraph model, as compared with 
the diffusion model. 

Finally, the telegraph equation proved to be a useful model in a wide range of transport 
problems (for a review, see Weiss 2002; Dunkel et ah 2007), and so the results of the present 
paper should have a broad applicability. 
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